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! Abstract 
^ . 

' Within a parton cascade we investigate the dependence of anisotropies in momentum space, 

• namely the elUptic flow V2 =< cos{2(j)) > and the =< cos{4:(j)) >, on both the finite shear 

^p^, viscosity rj and the freeze-out (f.o.) dynamics at the RHIC energy of 200 AGeV. In particular it 

is discussed the impact of the f.o. dynamics looking at two different procedures: switching-off the 
collisions when the energy density goes below a fixed value or reducing the cross section according 
to the increase in 77/5 from a QGP phase to a hadronic one. We address the relation between 
I I the scaling of 1^2 (pr) with the eccentricity €x and with the integrated elliptic flow. We show that 

^ ' the breaking of the V2{pt)/^x scaling is not coming mainly from the finite rj/s but from the f.o. 

. dynamics and that the V2{pt) is weakly dependent on the f.o. scheme. On the other hand the 

' 1^4 (pr) is found to be much more dependent on both the rj/s and the f.o. dynamics and hence is 

C*^ ' indicated to put better constraints on the properties of the QGP. A first semi-quantitative analysis 

y—i • show that both V2 and V4 (with the smooth f.o.) consistently indicate a plasma with Airrj/ s ~ 1 — 2. 

00 ' 

^ 1 Introduction 

> 

^ The first stage at the Relativistic Heavy Ion Collider (RHIC) has successfully shown the formation of 
hot and dense matter which behaves like a nearly perfect fluid. Such a conclusion mainly relies on the 
large value of the elliptic flow f 2 [H E] • However there is a growing evidence of a breakdown of the ideal 
hydrodynamical behavior especially in the intermediate region (1.5 < < 5 GeV) [21 HI [5] and in 
peripheral collisions that calls for a kinetic approach to study finite shear viscosity and non-equilibrium 
effects. A currently debated topic is in particular the effect of the conjectured minimal shear viscosity 
to entropy density ratio rj/s > l/Air [6J on collective flows. Several exploratory studies within viscous 
hydrodynamics [71 l8l iQl fTO] or cascade approaches [UlilSliTS] has shown a significant effect of the finite 
rj/s even for values around the lower bound, indicating a value for r)/s ~ 0.1 — 0.2. 

The analysis presented here is based on a parton cascade approach for massless particles and it 
is mainly focused on the intermediate pt region. The main idea is to keep the r]/s oi the medium 
constant during the collision dynamics by rearranging the parton cross section according to the local 
density and mean momentum values as suggested by D. Molnar in [Hj. We study the scaling of the 
V2{pt) with the eccentricity = {y"^ — x^) / (x^ -|- y^) as a function of the centrality of the collision for 
Au + Au collisions at ^/snn = 200 GeV. We show that if rj/s, in the range 1 < Anrj/s < 4, is kept 
constant down to the thermal freeze-out (e ~ 0.2 GeV/fm^) the scaling is not broken by the finite 
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shear viscosity in the whole px range investigated (up to 3.5 GeV). Hence the scahng is not a unique 
feature of ideal hydrodynamics [15j. On the other hand experimentally it has been observed a scaling of 
V2{pt) I {^2) [IS] with the centrality of the collision and the system size together with the broken scaling 
of V2{pt)I^x [HI- The question is what is the origin of the breaking once a value of 77/5 ~ 0.1 — 0.2 
is not found to be responsible for it. In Ref. [12] we pointed out that once a freeze-out condition is 
introduced at e/.o. ~ 0.7 GeV/fm^ a cascade approach at finite viscosity can account for the breaking of 
the scaling for V2{pt)/^x together with a persisting scaling for V2{pt)/ {^2), as experimentally observed. 
Here we present a further developments of such a study replacing the cut-off freeze-out based on a 
smooth increase of the r]/s in the region of the mixed quark-hadron phase (e ~ 0.3 — 1.7 GeV/fm^). 
We see that the V2{pt) is similar in the two freeze-out prescriptions confirming our first results |12j . 
Moreover we analyze for the first time the effect of finite 7]/s on the ^4 that is revealed to be much 
stronger than for the elliptic flow and much more sensitive to the way the system is described at energy 
densities e ~ 0.5 — 2 GeV/fm^. We finally emphasize that the shape of V2{pt) and V4^{pt) at Anrj/s ~ 2 
is consistent with to the one conjectured in coalescence models P|, hence a definitive evaluation of 77/5 
is entangled with the observation of quark number scaling in the same px range. Hadronization by 
coalescence plus fragmentation has to be self-consistently included in the next future. 



2 The parton cascade 

Our approach to study the effect of finite shear viscosity on the development of anisotropies in momen- 
tum space is based on a 3 + 1 dimensional Montecarlo cascade |12| for on-shell partons based on the 
stochastic interpretation of the transition rate thoroughly discussed in Ref. jl8j. Therefore the evolution 
of parton distribution function from initial conditions is governed by the Boltzmann equation 

P,d^fl= ///(/l'/2' - /l/2)|-Ml'2'^12p5'bl +P2 -Pi (1) 
21'2' 

where Jj = Jj cPpj/ [{211)^ 2Ej], M. denotes the transition matrix for the elastic processes and fj are 
the particle distribution functions. For the numerical implementation, we discretize the space into cells 
small respect to the system size and we use such cells to calculate all the local quantities. Several checks 
have been performed as in [18] to test the validity of the code and to choose a good discretization for 
the convergence of the results for momentum anisotropies. 

In kinetic theory under ultra-relativistic conditions the shear viscosity can be expressed as [19] 

ri=^p{p)\ (2) 

with p the parton density, A the mean free path and {p) the average momentum. Therefore considering 
that the entropy density for a massless gas is s = p(4 — p/T), fi being the chemical potential, we get: 

15 a,,p(4 - /i/T) 

where atr is the transport cross section, i.e. the sin^^ weighted cross section. We use a pQCD inspired 
cross section with the infrared singularity regularized by Debye thermal mass rriu [20] : 



da Qttq!? / 1 



D 



dt (t + ml)^\2 2t 
where s, t are the Mandelstam variables and niD = 0.7 GeV. 



(4) 
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Our approach is to artificially keep the rj/s of the medium constant during the dynamics of the 
collisions in a way similar to [HI [T3] , but evaluating locally in space and time the strength of the cross 
section atr{p{r),T) needed to keep the rj/s constant. From Eq. we see that assuming locally the 
thermal equilibrium this can be obtained evaluating in each a cell the cross section according to: 

_ 4 {p)a 1 

15p„(4-/.„/T)Vs 

with Attti/ s set in the range 1 — 4. We notice that a guideline on the temperature and time dependence of 
the cross section can be obtained considering the simple case of a free massless gas for which s = g^T^, 
and therefore neglecting fi in Eq. ([5]) one gets atr ~ for Anrj/s = 1. Furthermore a simple Bjorken 
expansion which means T ~ r^^/^ gives atr oc t"^^^ which is the approximate prescription adopted in 
[miTS]. In Fig{T]it is shown atr{T) evaluated locally in space in a cylinder of radius 3 fm as a function 
of time, we see on the left the approximate r^/^ and on the right the agreement with the estimated 
behavior. 




t [fm/c] T [GeV] 

Figure 1: Left: Time dependence of the cross section in a central region of rapidity {\y\ < 0.2) and 
transverse radius r < 3 fm. Right: Temperature dependence of the cross section in the central region 
compared with the dependence. 



3 Results 

Partons are initially distributed according to the Glauber model, hence also the eccentricity is similar 
to the one used in standard calculations. We also start our simulation like in hydrodynamics at a time 
t = 0.6 fm, assuming free- streaming evolution from t = to to t = 0.6/m. Partons with pt < Po =2 
GeV are distributed according to a thermalized spectrum, while for pt > Po we take the spectrum of 
non-quenched minijets as calculated in ^21j . 

A first important issue has been to investigate if the scaling behavior of the elliptic flow with the 
initial eccentricity and the system size typical of ideal hydrodynamics [IS] persists at finite t]/s close 
to the lower bound. The interest for such a study is triggered by the recent observation by PHENIX 
[T6] of a scaling of v^ipr)/ (^2) up to ~ 3 GeV, a region usually considered out of the range where 
hydrodynamics should work. In Ref. [12j it has been shown that both V2{pt)/^x and V2{pt)/{v2) scale 
with the impact parameter and system size if the fireball is let evolving down to vanishing energy density 
and a constant r]/s. This was found for both Au + Au and Cu + Cu at ^/s = 200 GeV. This indicates 
that the scaling V2{pt) / {^2), which is advocated as a signature of the hydrodynamical behavior [16], is a 
more general property that holds also at finite mean free path or shear viscosity at least for values close 
to the lower bound. Moreover the scaling is shown to persist also at higher pt (~ 3 GeV) where not 
only the scaling but also the saturation shape is correctly reproduced by the parton cascade approach. 
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We mention that our results are in qualitative agreement with predictions from viscous hydrodynamics 
[Hlin], however we are not aware of an explicit investigation of V2{pt)/^x and V2{pt)/ (^2) scaling within 
hydrodynamics. 

However this first result is in contrast with the experimental observation that V2{pt)/ {^2) scaling is 
observed while the V2{pt)/^x is not [17J. However it should be noticed that when the local energy 
density e ~ 1 GeV/fm^ the QGP phase stops and is likely that the V2 does not develops significantly 
during the hadronic phase as suggested also by recent analysis of experimental data [231 Sj- In Ref.[T2] 
for investigating the impact of a QGP freeze-out on the eccentricity scaling we used a cut-off freeze-out 
stopping the collisions among partons as the local energy density drops below e/.o. =0.7 GeV/fm^ an 
intermediate value in the range corresponding to a mixed quark-hadron phase 

3.1 Effect of QGP freeze-out on V2 and 

When the freeze-out condition is implemented a sizeable reduction for the elliptic flow is observed (see 
Fig. 2 in Ref. [12j). especially for the most peripheral collisions and at intermediate pt- Correspondingly 
the scaling of elliptic flow with the initial spatial eccentricity is broken [12]. In particular V2{pt)/^x 
varies of nearly 40% from b=3 fm to b=9 fm in the intermediate pt region (~3 GeV). The amount of 
such a spreading is consistent with the data reported by [T7j for the centrality selections — 10% and 
10 — 40%, with central collisions exhibiting a bigger elliptic flow to eccentricity ratio than the peripheral 
ones, see also Fig 121 On the other hand, the scaling of V2/ (^2) with the impact parameter is still present. 
We are therefore driven to the conclusion that the breaking of the V2{pt)/^ scaling traces back to the 
freeze out physics, which deserves a deeper investigation. 

This is partially due to the fact that {v2)/ex does not exactly scale with centrality, as pointed out 
by STAR and PHOBOS [31 [17] measurements and also qualitatively confirmed by our approach [12] . 
In hydrodynamics it has also recently been pointed out that the scaling with eccentricity for the (^2) is 
indeed slightly broken by the freeze-out [lOj . This is roughly speaking in line with our results. However 
there the effect is smaller because in hydrodynamics the viscosity is kept fixed at the minimum value 
also in the hadronic Dhase. 




Figure 2: Left: Behavior of the Atttj/s versus the energy density in the freeze-out procedures. The 
value for the hadrons gas is taken from [22]. Right: parton ^ (left panel) and -j^^ (right panel) in the 
central rapidity region {\y\ < 0.35) for Au-|-Au at -y/i = 200AGeV. Different symbols refer to cascade 
simulations at various impact parameters for Anri/s = 1. The results with the cut-off f.o. (dashed lines, 
open symbols) are compared with the one with the new f.o shown by solid lines and filled symbols. 

Our previous work drew the attention to the freeze-out (f.o.) but certainly its implementation as a 
sudden cut-off f.o. at the fixed e/.o. = 0.7 GeV/fm'^ is an oversimplified approach to the problem. Here 
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we present a first improvement that is realized changing the local t]/ s from the conjectured value in the 
QGP phase to the one estimated for an hadronic phase ^22j that is about eight times larger. The local 
ry/s starts to change when e < 1.7 GeV/fm^ and rises linearly with e reaching the values of Anrj/s = 8 
when e = 0.4 GeV/fm^ as shown in the left panel of Figl2] (left). Such a range is the one corresponding 
to the mixed phase of the standard equation of state (EOS) used in hydrodynamical calculations [1]. 
This procedure of course is more realistic respect to the sudden cut-off, however in Figl2] (right) we can 
see that the results are quite similar in the two f.o. schemes, i.e. with the cut-off at e/.o. (open symbols, 
dashed line) and for the new f.o. with the smooth increase of viscosity (filled symbols, solid line). 
The results of Ref. |12j and the relative discussion are therefore confirmed. Of course it is important 
to understand the physical motivation for the weak dependence of the V2{pt) on the f.o. procedure. 
Looking at the time scale one can realize that indeed most of the elliptic flow is formed earlier than 
the time scale at which the system is at e ~ e/.o. and hence the details of the interaction during such a 
phase are less relevant. In this perspective we think that this is also the origin of the moderate impact 
of different EOS on V2{pt) seen within hydrodynamic CZI 




Figure 3: Left: Differential elliptic flow for Au+Au collisions y/s^N =200 AGeV for two values of rj/s 
and the two f.o. schemes: the cut-off at e/.o. (dashed lines) and the new f.o. (solid line); squares are 
data from [16]. Right: v^lpr) for three values of rj/s in the new f.o. procedure (solid line) and the 
cut-off one (dashed line); data are from [26j. 



The study of azimuthal anisotropy can be extended to higher harmonics and in particular a finite 
fourth harmonic V4 =< cos{(j)) > [26j| has been measured at RHIC for the first time in heavy-ion 
collisions. A first study with AMPT model has shown that it is more sensible to the partonic dynamics 
|27j . We point out that it has a formation time that is delayed respect to the V2 one and therefore it 
becomes much more dependent on the dynamics at energy densities in the region of the mixed phase. 
This is clearly shown in Fig 13] where we can see that V2{pt) is essentially not affected by the f.o. scheme 
while V4{pt) drastically depends on it (compare the solid and dashed line both at Atttj/s =1). Moreover 
only the smooth f.o. procedure is able to reproduce a V4{pt) similar to the one observed in experiments 
(squares). Therefore we warn about the use of a cut-off f.o. as done also by other groups [11] even if 
the evaluation of the V2{pt) can still be considered quite reliable. The observation that v^^pt) is formed 
at energy density e ~ 0.5-2 GeV/fm^ let us envisage that and not V2 is a good probe of the EOS in 
the cross-over region and it would be interesting to examine it also within hydrodynamics. 

A last point we emphasize is the dependence of v^Ipt) on the value of rj/s. In FigHj there is a 
comparison between the relative change in V2{pt) (left) and V4,{pt) (right) when the Airrj/s is increased 
by a factor of two and four respect to the lower bound. It is seen that the v^Ipt) is more sensitive and 
therefore can in principle give a better constraint on the rj/s of the QGP. Moreover it seems, see FigJS] 
that also the v^Ipt) confirms that 7]/s is at most a couple of time the lower bound. 

Furthermore as for the elliptic flow there is a good agreement with data at low px while at interme- 
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Figure 4: Au+Au ^/sNN=200 AGeV and b=7fm. Left: relative enhancement of f2(pr) at Anrj/s =2 
(circles) and 4 (squares) respect to the one at 1. Right the same left panel but for the v^. 

diate px a room is open for the presence of a coalescence mechanism that would make the estimate of 
ri/s consistent with the observation of quark number scaling for both V2 and V4 [51 l^j . 



4 Conclusion 

We have investigated the dependence on the shear viscosity and the f.o. dynamics of the elliptic flow 
V2{pt) and of the v^^pt)- As a first result we find that the approximate scaling of V2{pt)/^x advocated 
as a signature of the perfect hydrodynamical behavior [TB] can still hold also at finite viscosity up 
to intermediate pt [I2]- However such a scaling versus centrality is present only if one makes the 
fireball evolve without any f.o. condition. We found in [T^ that the freeze-out condition for the 
partonic dynamics is essential to reproduce the sizeable breaking of the V2{pt)/^x scaling together with 
a persisting V2{pt)/ (^2) scaling. This indicates the relevance of the freeze-out of QGP dynamics. In a 
cascade approach it is frequently used a f.o. condition based on the switch-off of the collisions when 
the energy density e is below a fixed value e/.o.. Here we have also studied the effect of a smooth f.o. 
procedure that takes into account the change of rj/s going from the QGP to the hadronic matter. We 
find that both the f.o. modelling give quite similar V2{pt)/^x and V2{pt)/{v2), hence also the estimate 
of rj/s based on the V2 obtained with the "cut-off" at e/.o. can be considered quite stable. On the other 
hand we found a quite strong dependence of the v^^pt) on both the f.o. procedure and the rj/s. This 
is due to its formation scale that is close to the life-time of the QGP phase. We find that the e/.o. 
cut-off is not able at all to reproduce the shape and even the sign of the observed f4(pr) while the more 
realistic smooth f.o. produces a shape similar to the experimental data. A first estimate shows that 
also V4^{pt) indicates a plasma with Anrj/s ~ 1 — 2 as the analysis of f2(pr) does [3 El [9l [TTl [12] which 
is an encouraging consistency. 

As a last point we point out that at intermediate pt there are several evidences for hadronization 
via coalescence and it has been shown that due to a coalescence mechanism the parton V2 translates 
into a nearly doubled hadron V2 |5]. On the other hand the f.o. condition in a parton cascade approach 
generates a shape of W2,4(pr) that agrees with the data at low pt and leaves the room for a coalescence 
enhancement at intermediate pt- Therefore a definite evaluation ofrj/s from V2{pt) and/or V4^{pt) needs 
a further development of the parton cascade approach that includes self-consistently the coalescence and 
fragmentation process in order to account also for the baryon-meson quark number scaling. Finally a 
quantitative estimate should also include the effect of a non vanishing e — in the cross-over transition 
and is known to lower the sound velocity and hence the V2 [SI [T5j . 
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